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. |_H ] Abstract. A preliminary result from Markarian 421 observa- (Punch et al., 1992). Extensive measurements of the gamma- 
tions in the energy region above 10 TeV with the CANGAROO- ray flux in the TeV energy region have been carried out by 
^ , II 1 m telescope is presented. In January 200 1 , the HEGRA several groups using the ground-based imaging Cerenkov tele- 

■ - - ' group reported that Markarian 421 had become very active, scope (Krennrich et al., 1999; Aharonian et al., 1999; Piron 
with flux levels up to 4 times that of the Crab Nebula. As et al., 2001). The energy spectrum up to 10 TeV has been 
a result, we observed Mkn 421 during six nights from Jan- measured by observations at large zenith angles, up to 60 
uary 24th to February 1st, and four nights from March 1st degrees, and can be well described by a power law with a 
to 4th. Observations were carried out using the very large spectral index — 2.5^--^ —3.0. 

zenith angle technique (-70 degree) and the energy thresh- it is generally believed that TeV gamma-rays from exti-a- 

old is estimated from Monte Carlo simulations to be around galactic sources will suffer from flie absorption due to photon- 
10 TeV. We have detected gamma-ray emission in this energy photon collisions with inter-galactic infrared radiation (Gould 
range. & Schreder et al., 1967; Stecker, de Jager & Salamon, 1992). 

However, at present, no cutoff due to this absorption has been 

detected in the Mkn 421 spectrum. In order to determine 
meaningful constraints on this expected cutoff, it is important 
to extend measurements of the spectrum to energies above 
10 TeV. 

Markarian421 (Mkn421, z=0.031) was the first extra-galactic 

source from which TeV gamma-ray emission was detected I" *is paper, a preliminary result from Mkn 42 1 observa- 
tions with the CANGAROO-II 10 m telescope is presented. 



1 Introduction 



Correspondence to: K.Okumura (okumura@icrr.u- We observed this target at very large zenith angles of about 
tokyo.ac.jp) 70 degree and searched for gamma-rays above 10 TeV. A 
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Table 1. Date, modified Julian day and effective observation hours 
for MarJcarian 421 observations. 
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Monte Carlo study of this observation mode is also described. 



2 Observation of Mkn 421 



The CANGAROO-II telescope is located near Woomera, South 
Austtaha (136°47'E, 31°06'S). A detailed explanation of the 
CANGAROO-11 telescope is given elsewhere (Tanimori et 
al., 1999; Mori et al, 1999, 2000, 2001; Kubo et al., 2000; 
Kawachietal.,2001). 

The telescope was pointed in the direction of Mkn421 for 
six nights from January 24th to February 1st, during which 
two nights were lost due to bad weather, and four nights from 
March 1st to 4th in 2001 . Information about these observing 
periods is summarized in Table 2. Mkn 421 was in a high 
state during these observation periods. In the first period, the 
TeV activity was monitored night by night by the HEGRA 
group (HEGRA , http://www-hegra.desy.de/mrk-421/). Their 
preliminary result shows that flux Ues in the range of 1~4 
times that of the Crab Nebula. Although there was no HEGRA 
measurement in the second period, Mkn 421 was stiU in a ^ ^^^^ 

high state according to the X-ray measurements of the RXTE/ASN^ . ^ , . 
group (RXTE/ASM , http://xte.mit.edu/ASMJc.html). 

Mkn 421 was observed at zenith angles of ~70 degrees as 
the elevation angle in culmination is 20.7 degrees at the site 
of CANGAROO-II telescope. Observations were scheduled 
when the source was above an elevation angle of 18 degrees. 
The observable time amounts to about two and a half hours 
per night. Data were taken on moonless and clear nights. An 
OFF source run, displaced in right ascension, followed or 
preceded the ON source run. An event trigger was registered 
when 3 ~ 5 individual pixels exceeded a threshold of about 
3 photoelectrons, depending on the night sky background. 



Fig. 1. TeV and X-ray activity measured by the 
HEGRA and RXTE/ASM groups, from (HEGRA , 
http://www-hegra.desy.de/mrk-421/) and (RXTE/ASM 
http://xte.imt.edu/ASMJc.html). The one-day averaged data 
is shown for lightcurve by RXTE/ASM. The observation period of 
the CANGAROO-II telescope is shown by the hatched region. 



are eliminated by the following procedure of image clean- 
ing. Pixels with pulse heights greater than ^^3.3 photoelec- 
tron, triggering within 30 nanoseconds of the timing center 
of the event and with more than three neighboring hit pixels 
are used. Final acceptance requires more than five hit pixels 
in each event. 

In order to reject data affected by cloud etc., data for this 
analysis were selected based on the following criteria: Data 
recorded under good sky conditions are used; The elevation 
angle of Mkn 421 was above 18.5 degrees; and data, which 
event rate after the image cleaning procedure is less than by 
2 r.m.s. from the averaged rate of the whole data, are elimi- 
nated. This cut does not exclude any reasonable short-term 



As a result, observations totalhng 14 hours for ON source 
and 13 hours for OFF source were selected. The average ele- 
vation angle of the selected data is 20.2 degree and its r.m.s. 
is 0.6 degree, thus data which is taken around culmination 
are selected. 

The selection of gamma-ray events is based on the stan- 
dard imaging analysis parameters; "Width", "Length", "Dis- 
tance" and "Alpha" (HiUas, 1982). In this analysis, we adopted 
the Likelihood method (Enomoto et al., 2001) for the selec- 
tion of gamma-ray events, instead of the conventional param- 
eterization method. The discrimination of gamma-ray events 
is based on the probability ratio Rprob> expressed as follows 



3 Analysis method 

The ADC and TDC information recorded for individual tubes 
are calibrated with the events triggered with a blue LED flasher 
before each observation. Flat fielding of the photomultiph- 
ers' gains and correction of TDC slewing effect are carried 
out using this data. Camera pixels triggered by sky noise 



p _ Probjj) 

Probi-f) + Prob{B.G.) 

where Prob{'-f) and Prob{B.G.) are calculated probabilities 
of image shape parameters for gamma-ray and background 
events, respectively. These are the products of individual 
probabilities for "Width", "Length" and "Distance", which 
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Fig. 2. Shower parameter distribution of gamma-ray Monte Carlo 
simulation (solid line) and background (dashed line). The "Alpha" 
distribution is obtained for selected events by the cut of image pa- 
rameters. 



are derived from the functionized distribution using gamma- 
ray simulations and and OFF source events. Energy depen- 
dence is taken into account in these probability functions. We 
adopted a cut of 0.56, as the signal-to-noise (S/N) ratio was 
maximized at this value in the study using simulated gamma- 
rays and background events. 



4 Monte Carlo simulation 
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Fig. 3. Effective area (left) and its product with E~^ ° spectrum 
(right). The energy threshold is estimated to be around 9.3 TeV, 
which is defined as the maximal point of this distribution. 



atively farther from the telescope due to the thicker atmo- 
sphere. This implies a higher energy threshold with a broader 
effective area (Sommers & Elbert, 1987; Tanimori et al., 1994, 
1998). The shower images focused on camera plane become 
smaller and the location of images moves closer to the target 
direction. The parameter values become smaller, compared 
with those of the observation near zenith, due to these ef- 
fects. The resolution of the pointing angle, "Alpha", also de- 
creases, but is still usable. Naturally the signal-to-noise ratio 
is smaller than for higher elevation observations, but not so 
small that it is impossible to discriminate against cosmic-ray 
protons. 

Fig 3 shows the effective area and relative event rate for 
gamma-rays, which is the product of the effective area times 
the assumed spectrum of E~^-°. The energy threshold is de- 
fined where the event rate becomes maximal and estimated 
to be 9.3 TeV. The averaged effective area above this energy 
threshold is calculated to be 5 x 10^ cm-^. It increases to 10^° 
cm^ as the gamma-ray energy increases. 



We have carried out a Monte Carlo simulation of the de- 
velopment of electro-magnetic shower cascade, generation 
and attenuation of Cerenkov photons, reflection of photons 
at the mirror surface and simulation of electrical hardware. 
GEANT simulation code (CERN, 1994) is used for the frame- 
work of particle tracking. Photon attenuation in the atmo- 
sphere, mirror reflectivity, and quantum efficiency of pho- 
tomultipliers are included as a function of Cerenkov photon 
wavelength. The measured point spread function and mir- 
ror reflectivity is used. Performance of ADC and TDC elec- 
tronics is also simulated, including accidental hits due to sky 
noise. 

Gamma-ray events are generated in the energy range of 
3'--^50 TeV with the assumed spectrum of E^'^-'^. The eleva- 
tion angle of each event is adjusted according to the distribu- 
tion of selected data so as to reduce the difference in energy 
scale between real data and simulation, since the energy scale 
is very sensitive for large zenith angle observations. 

Fig 2 shows the parameter distributions obtained from gamma- 
ray simulations and OFF source events. For observations 
at large zenith angle, the shower development point is rel- 



5 Result and discussion 

Fig 4 shows the distribution of the image orientation angle 
"Alpha". The normahzation factor between ON and OFF 
observations is obtained from the number of events in the 
"Alpha" distributions between 40° and 90°. The gamma-ray 
signal region is alpha < 20 degrees, and the excess events 
after subtraction of the background is 221 ± 39 events, cor- 
responding to a significance of 5.6 a. The broader alpha 
acceptance is due to the shrinkage of the shower images as 
predicted from gamma-ray simulations at very large zenith 
angles. 

This preliminary result implies the detection of gamma- 
rays from Mkn 421 in the energy region of 10~50 TeV. A 
refined analysis will be important for studies of the infrared 
radiation absorption, as discussed in section 1 . Furthermore, 
we have demonstrated that there is sensitivity for the higher 
energy region above 10 TeV from the study of Monte Carlo 
simulation. Further analysis may enable constraints to be 
placed on the intergalactic infra-red background. 
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Fig. 4. Image orientation angle "Alpha" 
distribution of Markarian 421 observed 
by CANGAROO-II 10 meter telescope. 
Dots with error bars are ON source data 
and the solid line is for the OFF source 
data. The background-subtracted distri- 
bution is shown in the inserted graph. 



6 Conclusion 

We have observed Markarian 421 with the CANGAROO-II 
10 m telescope for 14 hours (ON source data after cuts) while 
the source was in a high state at both X-ray and TeV energies 
around the end of January and the beginning of March, 2001 . 
The observations were carried out at very large zenith angles, 
around 70 degrees. The viability of gamma-ray observa- 
tions are verified with Monte Carlo simulations. The energy 
threshold and effective area are estimated to be 9.3 TeV and 
5 X 10^ cm^, respectively. An excess of 221 ± 39 events was 
found using Likelihood method. This corresponds to a sig- 
nificance of 5.6 cr and the preliminary detection of gamma- 
rays in the energy range of 10~50 TeV. 
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